Autoimmunity Against Surfactant Protein B Is Associated with Pneumonitis During Checkpoint Blockade.

RATIONALE: Immune checkpoint inhibitor-related pneumonitis is a serious autoimmune event affecting up to 20% of patients with non-small cell lung cancer, yet the factors underpinning its development in some patients and not others are poorly understood. OBJECTIVES: To investigate the role of autoantibodies and autoreactive T cells against surfactant-related proteins in the development of pneumonitis. METHODS: The study cohort consisted of non-small cell lung cancer patients who gave blood samples before and during immune checkpoint inhibitor treatment. Serum was used for proteomics analyses and to detect autoantibodies present during pneumonitis. T cell stimulation assays and single-cell RNA sequencing were performed to investigate the specificity and functionality of peripheral autoreactive T cells. The findings were confirmed in a validation cohort comprising patients with non-small cell lung cancer and patients with melanoma. MEASUREMENTS AND MAIN RESULTS: Across both cohorts, patients who developed pneumonitis had higher pre-treatment levels of immunoglobulin G autoantibodies targeting surfactant protein-B. At the onset of pneumonitis, these patients also exhibited higher frequencies of CD4+ interferon-gamma-positive surfactant protein B-specific T cells, and expanding T cell clonotypes recognizing this protein, accompanied by a pro-inflammatory serum proteomic profile. CONCLUSIONS: Our data suggest that the co-occurrence of surfactant protein-B-specific immunoglobulin G autoantibodies and CD4+ T cells is associated with the development of pneumonitis during ICI therapy. Pre-treatment levels of these antibodies may represent a potential biomarker for elevated risk of developing pneumonitis and on-treatment levels may provide a diagnostic aid.

An immunogenetic basis for lung cancer risk.

Cancer risk is influenced by inherited mutations, DNA replication errors, and environmental factors. However, the influence of genetic variation in immunosurveillance on cancer risk is not well understood. Leveraging population-level data from the UK Biobank and FinnGen, we show that heterozygosity at the human leukocyte antigen (HLA)-II loci is associated with reduced lung cancer risk in smokers. Fine-mapping implicated amino acid heterozygosity in the HLA-II peptide binding groove in reduced lung cancer risk, and single-cell analyses showed that smoking drives enrichment of proinflammatory lung macrophages and HLA-II+ epithelial cells. In lung cancer, widespread loss of HLA-II heterozygosity (LOH) favored loss of alleles with larger neopeptide repertoires. Thus, our findings nominate genetic variation in immunosurveillance as a critical risk factor for lung cancer.

Template-specific optimization of NGS genotyping pipelines reveals allele-specific variation in MHC gene expression.

Using high-throughput sequencing for precise genotyping of multi-locus gene families, such as the major histocompatibility complex (MHC), remains challenging, due to the complexity of the data and difficulties in distinguishing genuine from erroneous variants. Several dedicated genotyping pipelines for data from high-throughput sequencing, such as next-generation sequencing (NGS), have been developed to tackle the ensuing risk of artificially inflated diversity. Here, we thoroughly assess three such multi-locus genotyping pipelines for NGS data, the DOC method, AmpliSAS and ACACIA, using MHC class IIß data sets of three-spined stickleback gDNA, cDNA and "artificial" plasmid samples with known allelic diversity. We show that genotyping of gDNA and plasmid samples at optimal pipeline parameters was highly accurate and reproducible across methods. However, for cDNA data, the gDNA-optimal parameter configuration yielded decreased overall genotyping precision and consistency between pipelines. Further adjustments of key clustering parameters were required tο account for higher error rates and larger variation in sequencing depth per allele, highlighting the importance of template-specific pipeline optimization for reliable genotyping of multi-locus gene families. Through accurate paired gDNA-cDNA typing and MHC-II haplotype inference, we show that MHC-II allele-specific expression levels correlate negatively with allele number across haplotypes. Lastly, sibship-assisted cDNA-typing of MHC-I revealed novel variants linked in haplotype blocks, and a higher-than-previously-reported individual MHC-I allelic diversity. In conclusion, we provide novel genotyping protocols for the three-spined stickleback MHC-I and -II genes, and evaluate the performance of popular NGS-genotyping pipelines. We also show that fine-tuned genotyping of paired gDNA-cDNA samples facilitates amplification bias-corrected MHC allele expression analysis.

Disruption of a Plasmodium falciparum patatin-like phospholipase delays male gametocyte exflagellation.

An essential process in transmission of the malaria parasite to the Anopheles vector is the conversion of mature gametocytes into gametes within the mosquito gut, where they egress from the red blood cell (RBC). During egress, male gametocytes undergo exflagellation, leading to the formation of eight haploid motile microgametes, while female gametes retain their spherical shape. Gametocyte egress depends on sequential disruption of the parasitophorous vacuole membrane and the host cell membrane. In other life cycle stages of the malaria parasite, phospholipases have been implicated in membrane disruption processes during egress, however their importance for gametocyte egress is relatively unknown. Here, we performed comprehensive functional analyses of six putative phospholipases for their role during development and egress of Plasmodium falciparum gametocytes. We localize two of them, the prodrug activation and resistance esterase (PF3D7_0709700) and the lysophospholipase 1 (PF3D7_1476700), to the parasite plasma membrane. Subsequently, we show that disruption of most of the studied phospholipase genes does neither affect gametocyte development nor egress. The exception is the putative patatin-like phospholipase 3 (PF3D7_0924000), whose gene deletion leads to a delay in male gametocyte exflagellation, indicating an important, albeit not essential, role of this enzyme in male gametogenesis.

Macrophage susceptibility to infection by Ghanaian Mycobacterium tuberculosis complex lineages 4 and 5 varies with self-reported ethnicity.

Background: The epidemiology of Mycobacterium tuberculosis complex (MTBC) lineage 5 (L5) infections in Ghana revealed a significantly increased prevalence in Ewes compared to other self-reported ethnic groups. In that context, we sought to investigate the early phase of tuberculosis (TB) infection using ex vivo infection of macrophages derived from the blood of Ewe and Akan ethnic group volunteers with MTBC L4 and L5 strains. Methods: The study participants consisted of 16 controls, among which self-reported Akan and Ewe ethnicity was equally represented, as well as 20 cured TB cases consisting of 11 Akans and 9 Ewes. Peripheral blood mononuclear cells were isolated from both healthy controls and cured TB cases. CD14+ monocytes were isolated and differentiated into monocyte-derived macrophages (MDMs) before infection with L4 or L5 endemic strains. The bacterial load was assessed after 2 hours (uptake) as well as 3 and 7 days post-infection. Results: We observed a higher capacity of MDMs from Ewes to phagocytose L4 strains (p < 0.001), translating into a higher bacillary load on day 7 (p < 0.001) compared to L5, despite the higher replication rate of L5 in Ewe MDMs (fold change: 1.4 vs. 1.2, p = 0.03) among the controls. On the contrary, within macrophages from Akans, we observed a significantly higher phagocytic uptake of L5 (p < 0.001) compared to L4, also translating into a higher load on day 7 (p = 0.04). However, the replication rate of L4 in Akan MDMs was higher than that of L5 (fold change: L4 = 1.2, L4 = 1.1, p = 0.04). Although there was no significant difference in the uptake of L4 and L5 among cured TB cases, there was a higher bacterial load of both L4 (p = 0.02) and L5 (p = 0.02) on day 7 in Ewe MDMs. Conclusion: Our results suggest that host ethnicity (driven by host genetic diversity), MTBC genetic diversity, and individual TB infection history are all acting together to modulate the outcome of macrophage infections by MTBC.

Tutorial: a statistical genetics guide to identifying HLA alleles driving complex disease.

The human leukocyte antigen (HLA) locus is associated with more complex diseases than any other locus in the human genome. In many diseases, HLA explains more heritability than all other known loci combined. In silico HLA imputation methods enable rapid and accurate estimation of HLA alleles in the millions of individuals that are already genotyped on microarrays. HLA imputation has been used to define causal variation in autoimmune diseases, such as type I diabetes, and in human immunodeficiency virus infection control. However, there are few guidelines on performing HLA imputation, association testing, and fine mapping. Here, we present a comprehensive tutorial to impute HLA alleles from genotype data. We provide detailed guidance on performing standard quality control measures for input genotyping data and describe options to impute HLA alleles and amino acids either locally or using the web-based Michigan Imputation Server, which hosts a multi-ancestry HLA imputation reference panel. We also offer best practice recommendations to conduct association tests to define the alleles, amino acids, and haplotypes that affect human traits. Along with the pipeline, we provide a step-by-step online guide with scripts and available software ( https://github.com/immunogenomics/HLA_analyses_tutorial ). This tutorial will be broadly applicable to large-scale genotyping data and will contribute to defining the role of HLA in human diseases across global populations.

Origin matters: Using a local reference genome improves measures in population genomics.

Genome sequencing enables answering fundamental questions about the genetic basis of adaptation, population structure and epigenetic mechanisms. Yet, we usually need a suitable reference genome for mapping population-level resequencing data. In some model systems, multiple reference genomes are available, giving the challenging task of determining which reference genome best suits the data. Here, we compared the use of two different reference genomes for the three-spined stickleback (Gasterosteus aculeatus), one novel genome derived from a European gynogenetic individual and the published reference genome of a North American individual. Specifically, we investigated the impact of using a local reference versus one generated from a distinct lineage on several common population genomics analyses. Through mapping genome resequencing data of 60 sticklebacks from across Europe and North America, we demonstrate that genetic distance among samples and the reference genomes impacts downstream analyses. Using a local reference genome increased mapping efficiency and genotyping accuracy, effectively retaining more and better data. Despite comparable distributions of the metrics generated across the genome using SNP data (i.e. π, Tajima's D and FST ), window-based statistics using different references resulted in different outlier genes and enriched gene functions. A marker-based analysis of DNA methylation distributions had a comparably high overlap in outlier genes and functions, yet with distinct differences depending on the reference genome. Overall, our results highlight how using a local reference genome decreases reference bias to increase confidence in downstream analyses of the data. Such results have significant implications in all reference-genome-based population genomic analyses.

Understanding the evolution of immune genes in jawed vertebrates.

Driven by co-evolution with pathogens, host immunity continuously adapts to optimize defence against pathogens within a given environment. Recent advances in genetics, genomics and transcriptomics have enabled a more detailed investigation into how immunogenetic variation shapes the diversity of immune responses seen across domestic and wild animal species. However, a deeper understanding of the diverse molecular mechanisms that shape immunity within and among species is still needed to gain insight into-and generate evolutionary hypotheses on-the ultimate drivers of immunological differences. Here, we discuss current advances in our understanding of molecular evolution underpinning jawed vertebrate immunity. First, we introduce the immunome concept, a framework for characterizing genes involved in immune defence from a comparative perspective, then we outline how immune genes of interest can be identified. Second, we focus on how different selection modes are observed acting across groups of immune genes and propose hypotheses to explain these differences. We then provide an overview of the approaches used so far to study the evolutionary heterogeneity of immune genes on macro and microevolutionary scales. Finally, we discuss some of the current evidence as to how specific pathogens affect the evolution of different groups of immune genes. This review results from the collective discussion on the current key challenges in evolutionary immunology conducted at the ESEB 2021 Online Satellite Symposium: Molecular evolution of the vertebrate immune system, from the lab to natural populations.

Challenges and opportunities in analyzing and modeling peptide presentation by HLA-II proteins.

The human leukocyte antigen (HLA) proteins are an indispensable component of adaptive immunity because of their role in presenting self and foreign peptides to T cells. Further, many complex diseases are associated with genetic variation in the HLA region, implying an important role for specific HLA-presented peptides in the etiology of these diseases. Identifying the specific set of peptides presented by an individual's HLA proteins in vivo, as a whole being referred to as the immunopeptidome, has therefore gathered increasing attention for different reasons. For example, identifying neoepitopes for cancer immunotherapy, vaccine development against infectious pathogens, or elucidating the role of HLA in autoimmunity. Despite the tremendous progress made during the last decade in these areas, several questions remain unanswered. In this perspective, we highlight five remaining key challenges in the analysis of peptide presentation and T cell immunogenicity and discuss potential solutions to these problems. We believe that addressing these questions would not only improve our understanding of disease etiology but will also have a direct translational impact in terms of engineering better vaccines and in developing more potent immunotherapies.

A Microtubule-Associated Protein Is Essential for Malaria Parasite Transmission.

Mature gametocytes of Plasmodium falciparum display a banana (falciform) shape conferred by a complex array of subpellicular microtubules (SPMT) associated with the inner membrane complex (IMC). Microtubule-associated proteins (MAPs) define MT populations and modulate interaction with pellicular components. Several MAPs have been identified in Toxoplasma gondii, and homologues can be found in the genomes of Plasmodium species, but the function of these proteins for asexual and sexual development of malaria parasites is still unknown. Here, we identified a novel subpellicular MAP, termed SPM3, that is conserved within the genus Plasmodium, especially within the subgenus Laverania, but absent in other Apicomplexa. Conditional knockdown and targeted gene disruption of Pfspm3 in Plasmodium falciparum cause severe morphological defects during gametocytogenesis, leading to round, nonfalciform gametocytes with an aberrant SPMT pattern. In contrast, Pbspm3 knockout in Plasmodium berghei, a species with round gametocytes, caused no defect in gametocytogenesis, but sporozoites displayed an aberrant motility and a dramatic defect in invasion of salivary glands, leading to a decreased efficiency in transmission. Electron microscopy revealed a dissociation of the SPMT from the IMC in Pbspm3 knockout parasites, suggesting a function of SPM3 in anchoring MTs to the IMC. Overall, our results highlight SPM3 as a pellicular component with essential functions for malaria parasite transmission. IMPORTANCE A key structural feature driving the transition between different life cycle stages of the malaria parasite is the unique three-membrane pellicle, consisting of the parasite plasma membrane (PPM) and a double membrane structure underlying the PPM termed the inner membrane complex (IMC). Additionally, there are numerous linearly arranged intramembranous particles (IMPs) linked to the IMC, which likely link the IMC to the subpellicular microtubule cytoskeleton. Here, we identified, localized, and characterized a novel subpellicular microtubule-associated protein unique to the genus Plasmodium. The knockout of this protein in the human-pathogenic species P. falciparum resulted in malformed gametocytes and aberrant microtubules. We confirmed the microtubule association in the P. berghei rodent malaria homologue and show that its knockout results in a perturbed microtubule architecture, aberrant sporozoite motility, and decreased transmission efficiency.

Intronic primers reveal unexpectedly high major histocompatibility complex diversity in Antarctic fur seals.

The major histocompatibility complex (MHC) is a group of genes comprising one of the most important components of the vertebrate immune system. Consequently, there has been much interest in characterising MHC variation and its relationship with fitness in a variety of species. Due to the exceptional polymorphism of MHC genes, careful PCR primer design is crucial for capturing all of the allelic variation present in a given species. We therefore developed intronic primers to amplify the full-length 267 bp protein-coding sequence of the MHC class II DQB exon 2 in the Antarctic fur seal. We then characterised patterns of MHC variation among mother-offspring pairs from two breeding colonies and detected 19 alleles among 771 clone sequences from 56 individuals. The distribution of alleles within and among individuals was consistent with a single-copy, classical DQB locus showing Mendelian inheritance. Amino acid similarity at the MHC was significantly associated with genome-wide relatedness, but no relationship was found between MHC heterozygosity and genome-wide heterozygosity. Finally, allelic diversity was several times higher than reported by a previous study based on partial exon sequences. This difference appears to be related to allele-specific amplification bias, implying that primer design can strongly impact the inference of MHC diversity.

Individual HLA heterogeneity and its implications for cellular immune evasion in cancer and beyond.

Structural and functional variability of human leukocyte antigen (HLA) is the foundation for competent adaptive immune responses against pathogen and tumor antigens as it assures the breadth of the presented immune-peptidome, theoretically sustaining an efficient and diverse T cell response. This variability is presumably the result of the continuous selection by pathogens, which over the course of evolution shaped the adaptive immune system favoring the assortment of a hyper-polymorphic HLA system able to elaborate efficient immune responses. Any genetic alteration affecting this diversity may lead to pathological processes, perturbing antigen presentation capabilities, T-cell reactivity and, to some extent, natural killer cell functionality. A highly variable germline HLA genotype can convey immunogenetic protection against infections, be associated with tumor surveillance or influence response to anti-neoplastic treatments. In contrast, somatic aberrations of HLA loci, rearranging the original germline configuration, theoretically decreasing its variability, can facilitate mechanisms of immune escape that promote tumor growth and immune resistance. The purpose of the present review is to provide a unified and up-to-date overview of the pathophysiological consequences related to the perturbations of the genomic heterogeneity of HLA complexes and their impact on human diseases, with a special focus on cancer.

Presence of autoantibodies in serum does not impact the occurrence of immune checkpoint inhibitor-induced hepatitis in a prospective cohort of cancer patients.

PURPOSE: Immune checkpoint inhibitor (ICI)-induced hepatitis belongs to the frequently occurring immune-related adverse events (irAEs), particularly with the combination therapy involving ipilimumab and nivolumab. However, predisposing factors predicting the occurrence of ICI-induced hepatitis are barely known. We investigated the association of preexisting autoantibodies in the development of ICI-induced hepatitis in a prospective cohort of cancer patients. METHODS: Data from a prospective biomarker cohort comprising melanoma and non-small cell lung cancer (NSCLC) patients were used to analyze the incidence of ICI-induced hepatitis, putatively associated factors, and outcome. RESULTS: 40 patients with melanoma and 91 patients with NSCLC received ICI between July 2016 and May 2019. 11 patients developed ICI-induced hepatitis (8.4%). Prior to treatment, 45.5% of patients in the hepatitis cohort and 43.8% of the control cohort showed elevated titers of autoantibodies commonly associated with autoimmune liver diseases (p = 0.82). We found two nominally significant associations between the occurrence of ICI-induced hepatitis and HLA alleles associated with autoimmune liver diseases among NSCLC patients. Of note, significantly more patients with ICI-induced hepatitis developed additional irAEs in other organs (p = 0.0001). Neither overall nor progression-free survival was affected in the hepatitis group. CONCLUSION: We found nominally significant associations of ICI-induced hepatitis with two HLA alleles. ICI-induced hepatitis showed no correlation with liver-specific autoantibodies, but frequently co-occurred with irAEs affecting other organs. Unlike other irAEs, ICI-induced hepatitis is not associated with a better prognosis.

Improved prediction of immune checkpoint blockade efficacy across multiple cancer types.

Only a fraction of patients with cancer respond to immune checkpoint blockade (ICB) treatment, but current decision-making procedures have limited accuracy. In this study, we developed a machine learning model to predict ICB response by integrating genomic, molecular, demographic and clinical data from a comprehensively curated cohort (MSK-IMPACT) with 1,479 patients treated with ICB across 16 different cancer types. In a retrospective analysis, the model achieved high sensitivity and specificity in predicting clinical response to immunotherapy and predicted both overall survival and progression-free survival in the test data across different cancer types. Our model significantly outperformed predictions based on tumor mutational burden, which was recently approved by the U.S. Food and Drug Administration for this purpose1. Additionally, the model provides quantitative assessments of the model features that are most salient for the predictions. We anticipate that this approach will substantially improve clinical decision-making in immunotherapy and inform future interventions.

The similarity of class II HLA genotypes defines patterns of autoreactivity in idiopathic bone marrow failure disorders.

Idiopathic aplastic anemia (IAA) is a rare autoimmune bone marrow failure (BMF) disorder initiated by a human leukocyte antigen (HLA)-restricted T-cell response to unknown antigens. As in other autoimmune disorders, the predilection for certain HLA profiles seems to represent an etiologic factor; however, the structure-function patterns involved in the self-presentation in this disease remain unclear. Herein, we analyzed the molecular landscape of HLA complexes of a cohort of 300 IAA patients and almost 3000 healthy and disease controls by deeply dissecting their genotypic configurations, functional divergence, self-antigen binding capabilities, and T-cell receptor (TCR) repertoire specificities. Specifically, analysis of the evolutionary divergence of HLA genotypes (HED) showed that IAA patients carried class II HLA molecules whose antigen-binding sites were characterized by a high level of structural homology, only partially explained by specific risk allele profiles. This pattern implies reduced HLA binding capabilities, confirmed by binding analysis of hematopoietic stem cell (HSC)-derived self-peptides. IAA phenotype was associated with the enrichment in a few amino acids at specific positions within the peptide-binding groove of DRB1 molecules, affecting the interface HLA-antigen-TCR ß and potentially constituting the basis of T-cell dysfunction and autoreactivity. When analyzing associations with clinical outcomes, low HED was associated with risk of malignant progression and worse survival, underlying reduced tumor surveillance in clearing potential neoantigens derived from mechanisms of clonal hematopoiesis. Our data shed light on the immunogenetic risk associated with IAA etiology and clonal evolution and on general pathophysiological mechanisms potentially involved in other autoimmune disorders.

Ancient DNA Study in Medieval Europeans Shows an Association Between HLA-DRB1*03 and Paratyphoid Fever.

Outbreaks of infectious diseases repeatedly affected medieval Europe, leaving behind a large number of dead often inhumed in mass graves. Human remains interred in two burial pits from 14th century CE Germany exhibited molecular evidence of Salmonella enterica Paratyphi C (S. Paratyphi C) infection. The pathogen is responsible for paratyphoid fever, which was likely the cause of death for the buried individuals. This finding presented the unique opportunity to conduct a paratyphoid fever association study in a European population. We focused on HLA-DRB1*03:01 that is a known risk allele for enteric fever in present-day South Asians. We generated HLA profiles for 29 medieval S. Paratyphi C cases and 24 contemporaneous controls and compared these to a modern German population. The frequency of the risk allele was higher in the medieval cases (29.6%) compared to the contemporaneous controls (13%; p = 0.189), albeit not significantly so, possibly because of small sample sizes. Indeed, in comparison with the modern controls (n = 39,689; 10.2%; p = 0.005) the frequency difference became statistically significant. This comparison also suggested a slight decrease in the allele's prevalence between the medieval and modern controls. Up to now, this is the first study on the genetic predisposition to Salmonella infection in Europeans and the only association analysis on paratyphoid fever C. Functional investigation using computational binding prediction between HLA variants and S. Paratyphi and S. Typhi peptides supported a reduced recognition capacity of bacterial proteins by DRB1*03:01 relative to other common DRB1 variants. This pattern could potentially explain the disease association. Our results suggest a slightly reduced predisposition to paratyphoid fever in modern Europeans. The causative allele, however, is still common today, which can be explained by a trade-off, as DRB1*03:01 is protective against infectious respiratory diseases such as severe respiratory syndrome (SARS). It is thus possible that the allele also provided resistance to corona-like viruses in the past.

Immunopeptidomics toolkit library (IPTK): a python-based modular toolbox for analyzing immunopeptidomics data.

BACKGROUND: The human leukocyte antigen (HLA) proteins play a fundamental role in the adaptive immune system as they present peptides to T cells. Mass-spectrometry-based immunopeptidomics is a promising and powerful tool for characterizing the immunopeptidomic landscape of HLA proteins, that is the peptides presented on HLA proteins. Despite the growing interest in the technology, and the recent rise of immunopeptidomics-specific identification pipelines, there is still a gap in data-analysis and software tools that are specialized in analyzing and visualizing immunopeptidomics data. RESULTS: We present the IPTK library which is an open-source Python-based library for analyzing, visualizing, comparing, and integrating different omics layers with the identified peptides for an in-depth characterization of the immunopeptidome. Using different datasets, we illustrate the ability of the library to enrich the result of the identified peptidomes. Also, we demonstrate the utility of the library in developing other software and tools by developing an easy-to-use dashboard that can be used for the interactive analysis of the results. CONCLUSION: IPTK provides a modular and extendable framework for analyzing and integrating immunopeptidomes with different omics layers. The library is deployed into PyPI at https://pypi.org/project/IPTKL/ and into Bioconda at https://anaconda.org/bioconda/iptkl , while the source code of the library and the dashboard, along with the online tutorials are available at https://github.com/ikmb/iptoolkit .

Unique Pathogen Peptidomes Facilitate Pathogen-Specific Selection and Specialization of MHC Alleles.

A key component of pathogen-specific adaptive immunity in vertebrates is the presentation of pathogen-derived antigenic peptides by major histocompatibility complex (MHC) molecules. The excessive polymorphism observed at MHC genes is widely presumed to result from the need to recognize diverse pathogens, a process called pathogen-driven balancing selection. This process assumes that pathogens differ in their peptidomes-the pool of short peptides derived from the pathogen's proteome-so that different pathogens select for different MHC variants with distinct peptide-binding properties. Here, we tested this assumption in a comprehensive data set of 51.9 Mio peptides, derived from the peptidomes of 36 representative human pathogens. Strikingly, we found that 39.7% of the 630 pairwise comparisons among pathogens yielded not a single shared peptide and only 1.8% of pathogen pairs shared more than 1% of their peptides. Indeed, 98.8% of all peptides were unique to a single pathogen species. Using computational binding prediction to characterize the binding specificities of 321 common human MHC class-I variants, we investigated quantitative differences among MHC variants with regard to binding peptides from distinct pathogens. Our analysis showed signatures of specialization toward specific pathogens especially by MHC variants with narrow peptide-binding repertoires. This supports the hypothesis that such fastidious MHC variants might be maintained in the population because they provide an advantage against particular pathogens. Overall, our results establish a key selection factor for the excessive allelic diversity at MHC genes observed in natural populations and illuminate the evolution of variable peptide-binding repertoires among MHC variants.

Transethnic analysis of the human leukocyte antigen region for ulcerative colitis reveals not only shared but also ethnicity-specific disease associations.

Inflammatory bowel disease (IBD) is a chronic inflammatory disease of the gut. Genetic association studies have identified the highly variable human leukocyte antigen (HLA) region as the strongest susceptibility locus for IBD and specifically DRB1*01:03 as a determining factor for ulcerative colitis (UC). However, for most of the association signal such as delineation could not be made because of tight structures of linkage disequilibrium within the HLA. The aim of this study was therefore to further characterize the HLA signal using a transethnic approach. We performed a comprehensive fine mapping of single HLA alleles in UC in a cohort of 9272 individuals with African American, East Asian, Puerto Rican, Indian and Iranian descent and 40 691 previously analyzed Caucasians, additionally analyzing whole HLA haplotypes. We computationally characterized the binding of associated HLA alleles to human self-peptides and analyzed the physicochemical properties of the HLA proteins and predicted self-peptidomes. Highlighting alleles of the HLA-DRB1*15 group and their correlated HLA-DQ-DR haplotypes, we not only identified consistent associations (regarding effects directions/magnitudes) across different ethnicities but also identified population-specific signals (regarding differences in allele frequencies). We observed that DRB1*01:03 is mostly present in individuals of Western European descent and hardly present in non-Caucasian individuals. We found peptides predicted to bind to risk HLA alleles to be rich in positively charged amino acids. We conclude that the HLA plays an important role for UC susceptibility across different ethnicities. This research further implicates specific features of peptides that are predicted to bind risk and protective HLA proteins.